BAR4R-PUB- 12/020 
SLAC-PUB-15134 



P. 
B 

C 



M 
C 



Evidence of B + — > t + u decays with hadronic B tags 

J. P. Lees, 1 V. Poireau, 1 V. Tisserand, 1 J. Garra Tico, 2 E. Grauges, 2 A. Palano ob , 3 G. Eigen, 4 B. Stugu, 4 
D. N. Brown, 5 L. T. Kerth, 5 Yu. G. Kolomensky, 5 G. Lynch, 5 H. Koch, 6 T. Schroeder, 6 D. J. Asgeirsson, 7 

C. Hearty, 7 T. S. Mattison, 7 J. A. McKenna, 7 R. Y. So, 7 A. Khan, 8 V. E. Blinov, 9 A. R. Buzykaev, 9 
V. P. Druzhinin, 9 V. B. Golubev, 9 E. A. Kravchenko, 9 A. P. Onuchin, 9 S. I. Serednyakov, 9 Yu. I. Skovpen, 9 
E. P. Solodov, 9 K. Yu. Todyshev, 9 A. N. Yushkov, 9 M. Bondioli, 10 D. Kirkby, 10 A. J. Lankford, 10 M. Mandelkern, 10 
H. Atmacan, 11 J. W. Gary, 11 F. Liu, 11 O. Long, 11 G. M. Vitug, 11 C. Campagnari, 12 T. M. Hong, 12 D. Kovalskyi, 12 
J. D. Richman, 12 C. A. West, 12 A. M. Eisner, 13 J. Kroseberg, 13 W. S. Lockman, 13 A. J. Martinez, 13 B. A. Schumm, 13 
A. Seiden, 13 D. S. Chao, 14 C. H. Cheng, 14 B. Echenard, 14 K. T. Flood, 14 D. G. Hitlin, 14 P. Ongmongkolkul, 14 

F. C. Porter, 14 A. Y. Rakitin, 14 R. Andreassen, 15 Z. Huard, 15 B. T. Meadows, 15 M. D. Sokoloff, 15 L. Sun, 15 
C. Bloom, 16 W. T. Ford, 16 A. Gaz, 16 U. Nauenberg, 16 J. G. Smith, 16 S. R. Wagner, 16 R. Ayad, 17 'Q W. H. Toki, 17 
Spaan, 18 K. R. Schubert, 19 R. Schwierz, 19 D. Bernard, 20 M. Verderi, 20 P. J. Clark, 21 S. Playfer, 21 D. Bettoni a , 22 
I Bozzi a , 22 R. Calabrese ab , 22 G. Cibinetto ab , 22 E. Fioravanti ab , 22 I. Garzia ab , 22 E. Luppi ab , 22 M. Munerato ab , 22 

L. Piemontese a , 22 V. Santoro a , 22 R. Baldini-Ferroli, 23 A. Calcaterra, 23 R. de Sangro, 23 G. Finocchiaro, 23 
Patteri, 23 I. M. Peruzzi, 23 'Q M. Piccolo, 23 M. Rama, 23 A. Zallo, 23 R. Contri ab , 24 E. Guido ab , 24 M. Lo Vetere ab , 24 
M. R. Monge ab , 24 S. Passaggio a , 24 C. Patrignani ab , 24 E. Robutti a , 24 B. Bhuyan, 25 V. Prasad, 25 C. L. Lee, 26 
. Morii, 26 A. J. Edwards, 27 A. Adametz, 28 U. Uwer, 28 H. M. Lacker, 29 T. Lueck, 29 P. D. Dauncey, 30 U. Mallik, 31 
C. Chen, 32 J. Cochran, 32 W. T. Meyer, 32 S. Prell, 32 A. E. Rubin, 32 A. V. Gritsan, 33 Z. J. Guo, 33 N. Arnaud, 34 

M. Davier, 34 D. Derkach, 34 G. Grosdidier, 34 F. Le Diberder, 34 A. M. Lutz, 34 B. Malaescu, 34 P. Roudeau, 34 
M. H. Schune, 34 A. Stocchi, 34 G. Wormser, 34 D. J. Lange, 35 D. M. Wright, 35 C. A. Chavez, 36 J. P. Coleman, 36 
J. R. Fry, 36 E. Gabathuler, 36 D. E. Hutchcroft, 36 D. J. Payne, 36 C. Touramanis, 36 A. J. Bevan, 37 F. Di Lodovico, 37 
R. Sacco, 37 M. Sigamani, 37 G. Cowan, 38 D. N. Brown, 39 C. L. Davis, 39 A. G. Denig, 40 M. Fritsch, 40 W. Gradl, 40 
K. Griessinger, 40 A. Hafner, 40 E. Prencipe, 40 R. J. Barlow, 41 '! G. Jackson, 41 G. D. Lafferty, 41 E. Behn, 42 
R. Cenci, 42 B. Hamilton, 42 A. Jawahery, 42 D. A. Roberts 42 C. Dallapiccola, 43 R. Cowan, 44 D. Dujmic, 44 

G. Sciolla, 44 R. Cheaib, 45 D. Lindemann, 45 P. M. Patel, 45 >| S. H. Robertson, 45 P. Biassoni ab , 46 N. Neri a , 46 
F. Palombo ab , 46 S. Stracka ab , 46 L. Cremaldi, 47 R. Godang, 47 <0 R. Kroeger, 47 P. Sonnek, 47 D. J. Summers, 47 
X. Nguyen, 48 M. Simard, 48 P. Taras, 48 G. De Nardo ab , 49 D. Monorchio ab , 49 G. Onorato ab , 49 C. Sciacca ab , 49 

M. Martinelli, 50 G. Raven, 50 C. P. Jessop, 51 J. M. LoSecco, 51 W. F. Wang, 51 K. Honscheid, 52 R. Kass, 52 
J. Brau, 53 R. Frey, 53 N. B. Sinev, 53 D. Strom, 53 E. Torrence, 53 E. Feltresi ab , 54 N. Gagliardi ab , 54 M. Margoni ab , 54 
M. Morandin , 54 M. Posocco a , 54 M. Rotondo a , 54 G. Simi a , 54 F. Simonetto ab , 54 R. Stroili ab , 54 S. Akar, 55 
E. Ben-Haim, 55 M. Bomben, 55 G. R. Bonneaud, 55 H. Briand, 55 G. Calderini, 55 J. Chauveau, 55 O. Hamon, 55 
Ph. Leruste, 55 G. Marchiori, 55 J. Ocariz, 55 S. Sitt, 55 M. Biasini ab , 56 E. Manoni ab , 56 S. Pacetti ab , 56 A. Rossi ab , 56 
C. Angelini Qb , 57 G. Batignani ab , 57 S. Bettarini ab , 57 M. Carpinelli ab , 57 >EI G. Casarosa ab , 57 A. Cervelli ab , 57 
F. Forti ab , 57 M. A. Giorgi ab , 57 A. Lusiani ac , 57 B. Oberhof ab , 57 E. Paoloni ab , 57 A. Perez a , 57 G. Rizzo ab , 57 
J. J. Walsh a , 57 D. Lopes Pegna, 58 J. Olsen, 58 A. J. S. Smith, 58 A. V. Telnov, 58 F. Anulli", 59 R. Faccini ab , 59 
F. Ferrarotto a , 59 F. Ferroni ab , 59 M. Gaspero ab , 59 L. Li Gioi Q 59 M. A. Mazzoni a , 59 G. Piredda a , 59 C. Biinger, 60 
O. Griinberg, 60 T. Hartmann, 60 T. Leddig, 60 H. Schroder, 60 ! C. Voss, 60 R. Waldi, 60 T. Adye, 61 E. O. Olaiya, 61 
F. F. Wilson, 61 S. Emery, 62 G. Hamel de Monchenault, 62 G. Vasseur, 62 Ch. Yeche, 62 D. Aston, 63 D. J. Bard, 63 
R. Bartoldus, 63 J. F. Benitez, 63 C. Cartaro, 63 M. R. Convery, 63 J. Dorfan, 63 G. P. Dubois-Felsmann, 63 
W. Dunwoodie, 63 M. Ebert, 63 R. C. Field, 63 M. Franco Sevilla, 63 B. G. Fulsom, 63 A. M. Gabareen, 63 
M. T. Graham, 63 P. Grenier, 63 C. Hast, 63 W. R. Lines, 63 M. H. Kelsey, 63 P. Kim, 63 M. L. Kocian, 63 
D. W. G. S. Leith, 63 P. Lewis, 63 B. Lindquist, 63 S. Luitz, 63 V. Luth, 63 H. L. Lynch, 63 D. B. MacFarlane, 63 
D. R. Muller, 63 H. Neal, 63 S. Nelson, 63 M. Perl, 63 T. Pulliam, 63 B. N. Ratcliff, 63 A. Roodman, 63 A. A. Salnikov, 63 
R. H. Schindler, 63 A. Snyder, 63 D. Su, 63 M. K. Sullivan, 63 J. Va'vra, 63 A. P. Wagner, 63 W. J. Wisniewski, 63 
M. Wittgen, 63 D. H. Wright, 63 H. W. Wulsin, 63 C. C. Young, 63 V. Ziegler, 63 W. Park, 64 M. V. Purohit, 64 
R. M. White, 64 J. R. Wilson, 64 A. Randle-Conde, 65 S. J. Sekula, 65 M. Bellis, 66 P. R. Burchat, 66 T. S. Miyashita, 66 
M. S. Alam, 67 J. A. Ernst, 67 R. Gorodeisky, 68 N. Guttman, 68 D. R. Peimer, 68 A. Soffer, 68 P. Lund, 69 
S. M. Spanier, 69 J. L. Ritchie, 70 A. M. Ruland, 70 R. F. Schwitters, 70 B. C. Wray, 70 J. M. Izen, 71 X. C. Lou, 71 
F. Bianchi Qb , 72 D. Gamba ab , 72 S. Zambito ab , 72 L. Lanceri ab , 73 L. Vitale ab , 73 F. Martinez- Vidal, 74 A. Oyanguren, 74 



2 



H. Ahmed, 75 J. Albert, 75 Sw. Banerjee, 75 F. U. Bernlochner, 75 H. H. F. Choi, 75 G. J. King, 75 R. Kowalewski, 75 
M. J. Lewczuk, 75 I. M. Nugent, 75 J. M. Roney, 75 R. J. Sobie, 75 N. Tasneem, 75 T. J. Gershon, 76 P. F. Harrison, 76 
T. E. Latham, 76 E. M. T. Puccio, 76 H. R. Band, 77 S. Dasu, 77 Y. Pan, 77 R. Prepost, 77 and S. L. Wu 77 

(The BABAR Collaboration) 

1 Laboratoire d'Annecy-le-Vieux de Physique des Particules (LAPP), 
Universite de Savoie, CNRS/IN2P3, F-74941 Annecy-Le-Vieux, France 
2 Universitat de Barcelona, Facultat de Fisica, Departament ECM, E-08028 Barcelona, Spain 
4 INFN Sezione di Bart ; Dipartimento di Fisica, Universitd di Barf , 1-70126 Bari, Italy 
4 University of Bergen, Institute of Physics, N-5007 Bergen, Norway 
5 Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 
f 'Ruhr Universitat Bochum, Institut fur Experimentalphysik 1, D-44^80 Bochum, Germany 
7 University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1 
8 Brunel University, Uxbridge, Middlesex UBS 3PH, United Kingdom 
g Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 
10 University of California at Irvine, Irvine, California 92697, USA 
11 University of California at Riverside, Riverside, California 92521, USA 
12 University of California at Santa Barbara, Santa Barbara, California 93106, USA 
13 University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 
14 California Institute of Technology, Pasadena, California 91125, USA 

15 University of Cincinnati, Cincinnati, Ohio 45221, USA 

16 University of Colorado, Boulder, Colorado 80309, USA 
"Colorado State University, Fort Collins, Colorado 80523, USA 

18 Technische Universitat Dortmund, Fakultdt Physik, D-44^21 Dortmund, Germany 
19 Technische Universitat Dresden, Institut fur Kern- und Teilchenphysik, D-01062 Dresden, Germany 
20 Laboratoire Leprince-Ringuet, Ecole Polytechnique, CNRS/IN2P3, F-91128 Palaiseau, France 
21 University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 
22 INFN Sezione di Ferrara a ; Dipartimento di Fisica, Universitd di Ferrara b , 1-44100 Ferrara, Italy 

2S INFN Laboratori Nazionali di Frascati, 1-00044 Frascati, Italy 
24 INFN Sezione di Genova a ; Dipartimento di Fisica, Universitd di Genova b , 1-16146 Genova, Italy 
25 Indian Institute of Technology Guwahati, Guwahati, Assam, 781 039, India 
26 Harvard University, Cambridge, Massachusetts 02138, USA 
27 Harvey Mudd College, Claremont, California 91711, USA 
28 Universitat Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 
29 Humboldt- Universitat zu Berlin, Institut fur Physik, Newtonstr. 15, D- 12489 Berlin, Germany 
30 Imperial College London, London, SW7 2AZ, United Kingdom 
31 University of Iowa, Iowa City, Iowa 52242, USA 
32 Iowa State University, Ames, Iowa 50011-3160, USA 
33 Johns Hopkins University, Baltimore, Maryland 21218, USA 
34 Laboratoire de I'Accelerateur Lineaire, IN2P3/CNRS et Universite Paris-Sud 11, 
Centre Scientifique d'Orsay, B. P. 34, F-91898 Orsay Cedex, France 
35 Lawrence Livermore National Laboratory, Livermore, California 94550, USA 
36 University of Liverpool, Liverpool L69 7ZE, United Kingdom 
37 Queen Mary, University of London, London, El 4NS, United Kingdom 
38 University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 
39 University of Louisville, Louisville, Kentucky 40292, USA 
41 Johannes Gutenberg- Universitat Mainz, Institut fur Kernphysik, D-55099 Mainz, Germany 
41 University of Manchester, Manchester M13 9PL, United Kingdom 

42 University of Maryland, College Park, Maryland 20742, USA 
43 University of Massachusetts, Amherst, Massachusetts 01003, USA 
44 Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, USA 
45 McGill University, Montreal, Quebec, Canada H3A 2T8 
46 INFN Sezione di Milano a ; Dipartimento di Fisica, Universitd di Milano b , 1-20133 Milano, Italy 
47 University of Mississippi, University, Mississippi 38677, USA 
48 Universite de Montreal, Physique des Particules, Montreal, Quebec, Canada H3C 3J7 
49 INFN Sezione di Napoli a ; Dipartimento di Scienze Fisiche, 
Universitd di Napoli Federico IF, 1-80126 Napoli, Italy 
50 NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands 

51 University of Notre Dame, Notre Dame, Indiana 46556, USA 
52 Ohio State University, Columbus, Ohio 43210, USA 
53 University of Oregon, Eugene, Oregon 97403, USA 
54 INFN Sezione di Padova a ; Dipartimento di Fisica, Universitd di Padova b , 1-35131 Padova, Italy 



3 



55 Laboratoire de Physique Nucleaire et de Hautes Energies, 
IN2P3/CNRS, Universite Pierre et Marie Curie-Paris6, 
Universite Denis Diderot- Paris!, F-75252 Paris, France 
56 INFN Sezione di Perugia? ; Dipartimento di Fisica, Universitd di Perugia b , 1-06100 Perugia, Italy 
57 INFN Sezione di Pisa '; Dipartimento di Fisica, 
Universitd di Pisa h ; Scuola Normale Superiore di Pisa c , 1-56127 Pisa, Italy 
58 Princeton University, Princeton, New Jersey 08544, USA 
INFN Sezione di Roma a ; Dipartimento di Fisica, 
Universitd di Roma La Sapienza , 1-00185 Roma, Italy 
60 Universitdt Rostock, D-18051 Rostock, Germany 
61 Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 
62 CEA, Irfu, SPP, Centre de Saclay, F-91191 Gif-sur-Yvette, France 
b:l SLAG National Accelerator Laboratory, Stanford, California 94309 USA 
6Jt University of South Carolina, Columbia, South Carolina 29208, USA 
65 Southern Methodist University, Dallas, Texas 75275, USA 
66 Stanford University, Stanford, California 94305-4060, USA 
67 State University of New York, Albany, New York 12222, USA 
68 Tel Aviv University, School of Physics and Astronomy, Tel Aviv, 69978, Israel 

69 University of Tennessee, Knoxville, Tennessee 37996, USA 

70 University of Texas at Austin, Austin, Texas 78712, USA 
University of Texas at Dallas, Richardson, Texas 75083, USA 

72 INFN Sezione di Torino 0, ; Dipartimento di Fisica Sperimentale, Universitd di Torino , 1-10125 Torino, Itab 
71 INFN Sezione di Trieste " ; Dipartimento di Fisica, Universitd di Trieste , T34127 Trieste, Italy 
7J, IFIC, Universitat de Valencia- CSIC, E- 46071 Valencia, Spain 
75 University of Victoria, Victoria, British Columbia, Canada V8W 3P6 
76 Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 
77 University of Wisconsin, Madison, Wisconsin 53706, USA 
(Dated: July 3, 2012) 

We present a search for the decay B + — > t+v using 467.8 x 10 6 BB pairs collected at the Y(4S) 
resonance with the BABAR detector at the SLAC PEP-II B-Factory. We select a sample of events 
with one completely reconstructed B~ in the hadronic decay mode (B~ — > D^°X~ and B~ — > 
J/ipX~). We examine the rest of the event to search for a B + — > t+u decay. We identify the t+ 
lepton in the following modes: t+ —¥ e+vV, t+ — > fj,+uV, t+ — > tx+v and t+ — > p+V. We find an 
excess of events with respect to the expected background, which excludes the null signal hypothesis 
at the level of 3.8cr (including systematic uncertainties) and corresponds to a branching fraction 
central value of B{B + t+u) = (1.83±g;g(stat.) ± 0.24(syst.)) x 10" 4 . 

PACS numbers: 13.20.-v, 13.25.Hw 



The study of the purely leptonic decay B+ — > t+u 
[l| is of particular interest to test the predictions of the 
Standard Model (SM) and to probe of new physics ef- 
fects. It is sensitive to the product of the B meson de- 
cay constant fs , and the absolute value of the Cabibbo- 
Kobayashi-Maskawa matrix element \V u b\ ,2j. In the SM 
the branching fraction is given by: 



B(B+ -> T + v) = 



G F m B m 2 
8V 



1 - 



JB\ 



b\ T B +, 



(1) 
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where Gf is the Fermi constant, m B and m T are, the 
B + meson and r lepton masses, respectively, and r B+ 
is the B + lifetime. Using the Lattice QCD calculation 
of f B = (189 ± 4) MeV 3, and the BABAR measure- 
ment of \V u b\ from charmless semileptonic B exclusive 
decays [4|, the predicted SM value of the brancing frac- 
tion is B SM (B+ -> t+u) = (0.62 ± 0.12) x 10" 4 . If 
we use the BABAR measurement of \V u b\ from inclusive 
charmless semileptonic B decays the SM prediction 
is B S m{B+ -> t+u) = (1.18 ± 0.16) x 10" 4 . 

The process is sensitive to possible extensions of 
the SM. For instance, in two-Higgs doublet mod- 
els (2HDM) Q and in minimal supersymmetric exten- 
sions [3| it can be mediated by a charged Higgs boson. 
A branching fraction measurement can, therefore, also 
be used to constrain the parameter space of new physics 
models. 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II storage ring. The sample 
corresponds to an integrated luminosity of 426 fb -1 at the 
T(4S) resonance. The sample contains (467.8±5.1) x 10 6 
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TABLE I: Published results for B + -> r + u from BABAR and 



Belle collaborations. 



Experiment Tag 


Branching Fraction ( x 10 4 ) 


BABAR 
BABAR 
Belle 
Belle 


hadronic [8] 
semileptonic [9] 
hadronic [10] 
semileptonic fill 


l-8±o l ± 0.4 ± 0.2 
1.7 ±0.8 ±0.2 

i 7Q +0.56+0.46 
i -' a -0.49-0.51 
-, r A +0.38 +0.29 
1 -°^-().'.ir-0.31 



BB decays (N B -p). The detector is described in detail 
elsewhere [liq . Charged particle trajectories are mea- 
sured in the tracking system composed of a five-layer 
double-sided silicon vertex tracker and a 40-layer drift 
chamber, operating in a 1.5 T solenoidal magnetic field. 
A Cherenkov detector is used for charged ir-K discrimi- 
nation, a Csl calorimeter for photon and electron identifi- 
cation, and the flux return of the solenoid, which consists 
of layers of iron interspersed with resistive plate chambers 
or limited streamer tubes, for muon and neutral hadron 
identification. 

We use a Monte Carlo (MC) simulation based on 
GEANT4 [l3| to estimate signal selection efficiencies and to 
study backgrounds. In MC simulated signal events, one 
B + meson decays as B + — > t + v and the other decays in 
any final state. The BB and continuum MC samples are 
equivalent to approximately 3 times and 1.5 times the 
data sample, respectively. Beam-related background and 
detector noise are sampled from data and overlaid on the 
simulated events. 

We reconstruct an exclusive decay of one of the B 
mesons in the event (which we refer to as the tag-i?) 
and examine the rest of the event for the experimen- 
tal signature of B + — > t + v. The tag-i? reconstruction 
can be performed by looking at both hadronic B decays 
and semileptonic B decays. Published results from both 
BABAR and Belle are summarized in Table |U 

We reconstruct the tag-i? candidate in the set of 

hadronic decays B~ — > M°X~, where M° denotes a 
D (*)o 

or a J/ip , and X denotes a system of hadrons 
with total charge —1 composed of rii7r ± , n-iK^ , 71371-°, 
n^K® where n\ + ri2 < 5, 7J2, and < 2. We recon- 
struct the D° as L>° i^7r + , K^tt+tt", K^tt + tt-tt + , 
K° s ir°, K q s tt + tt-, K Q q Tr+TT-n Q , K + K~, or tt+tt-. We 
reconstruct the D*° meson as D* a — > D a Tr°,D°j, and 
the J/tp meson via their decays J/ip — > e + e~,/i + /i~. 
Two kinematic variables are used to discriminate be- 
tween correctly reconstructed tag-i? candidates and mis- 
reconstructed events: the beam energy-substituted mass 
toes = \/ s /4 — p B , and the energy difference AE = 
Eb — V~s/2, where y/s is the total energy in the T(4S') 
center-of-mass system (CM) and ps and Eb respectively 
denote the momentum and the energy of the tag--B can- 
didate in the CM. The resolution on AE is measured to 
be a = 10 — 35 MeV, depending on the decay mode; 
we require \AE\ < 3cr^^. Events with a tag-i? can- 
didate arise from two possible classes with different Toes 



distributions. One class includes signal events with a cor- 
rectly reconstructed tag-i?, and background events from 
T(AS) —> B + B~ with a correctly reconstructed tag-i?. 
All these events are characterized by an toes distribution 
peaked at the nominal B mass (signal and peaking back- 
ground) . The other classes of events consist of continuum 
background, e + e~ — > qq (q = u, , s, ) and e + e~ — » t + t~, 

and combinatorial background, T(4S*) — > B°B° or B + B~ 
in which the tag-i? is misreconstructed. These events are 
characterized by a smooth toes distribution. 

If multiple tag-i? candidates are reconstructed in the 
event, we select the one with the lowest value of \ AE\. Af- 
ter the reconstruction of the tag-i?, we require the pres- 
ence of only one well-reconstructed track (signal track), 
with charge opposite to that of the tag-i?. The pu- 
rity V of each reconstructed tag-i? decay mode is esti- 
mated as the ratio of the number of peaking events with 
m ES > 5.27 GeV to the total number of events in the 
same range. The yield in data is determined by means 
of an extended unbinned maximum likelihood fit to the 
toes distribution, as shown in Fig. [TJ We use a phe- 
nomenologically motivated threshold function (ARGUS 
function [lj]) as probability density function (PDF) to 
describe the continuum and combinatorial background 
components in the fit, while for the correctly recon- 
structed tag-i? component we use a Gaussian distribu- 
tion plus an exponential tail for the PDF (Crystal Ball 
function) [l5| . We use only events with the tag-i? recon- 
structed in decay modes with V > 0.1. Combinatorial 
and continuum background distributions in any discrim- 
inating variable are estimated from a sideband in toes 
(5.209 GeV < to E s < 5.260 GeV) and are extrapolated 
into the signal region (toes > 5.270 GeV) using the re- 
sults of a fit to an ARGUS function. The peaking B + B~ 
background shape is determined from B + B~ MC, after 
subtraction of the combinatorial component to avoid dou- 
ble counting. 




m ES [GeV] 



FIG. 1: Fit to the mES distribution in data. Dots are data, 
the upper curve is the global fit result and the lower curve rep- 
resents the fitted combinatorial and continuum background. 
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The signal-side r lepton is reconstructed in four de- 
cay modes: r + — > e + vD, r + — > p^vv, r + — > tt + v 1 and 
t + — > totaling approximately 70% of all r decays. 
We separate the event sample into four categories using 
particle identification criteria applied to the signal track 
(e + , and 7r + ). The t + p + v sample is obtained by 
associating the signal track 7r + with a n° reconstructed 
from a pair of neutral clusters with an invariant mass 
between 115 MeV/c 2 and 155 MeV/c 2 . 

In order to remove the e + e~ — > t + t~ background, we 
impose r mode dependent requirements on the ratio be- 
tween the 2 nd and the 0*' 1 Fox- Wolfram moments R2 [16J 
calculated using all the tracks and neutral clusters of the 
event. This preserves 90% of the B + — > t + v signal. 

To reject continuum background, we use the absolute 
value of cos 9tb, the cosine of the angle in the CM frame 
between the thrust axis [17[ of the tag-B and the thrust 
axis of the remaining charged and neutral candidates 
in the event. For correctly reconstructed tag-B candi- 
dates the Icos^tbI distribution is expected to be uniform, 
while for jet-like e + e~ — > qq continuum events it peaks 
strongly at 1. In order to reject background from events 
with a correctly reconstructed tag-B, we study the distri- 
bution of several discriminating variables exploiting the 
different kinematics between the signal and background 
of the remaining reconstructed candidates. We use the 
missing momentum polar angle in the laboratory frame 
Pmiss = Pcm - PtagB - Ptrk - S„ e „ t Pij where pen is 
the total momentum of the beams, ptagB is the recon- 
structed momentum of the tag-B, and p tr k is the recon- 
structed track momentum, and the sum is extended on 
all the neutral candidates reconstructed in the calorime- 
ter not assigned to the tag-B. For the r + — > tt + v mode, 
we combine p* trk (where the star denotes the CM frame) 
and the cosine of the angle between p m %ss and the beam 
axis (cos Bmiss) m a likelihood ratio 

Lp = L s (p* trk ,cos6 miss ) 

{Lsipl rk , cos 6 

miss 

) + L B {p* trk ,cosO miss )y 

where the signal (S) and background (B) likeli- 
hoods have been obtained from the product of 
the PDFs of the two discriminating variables: 

Ls(P* trk iCOs0 miss ) = Ps(P$rk) P S (COS miss ) and 

L B{p* trk , cos miss ) = PB{Pt rk )PB (cos 9 miss ). Sim- 
ilarly, for the t + — > p + v mode we combine four 
discriminating variables in the likelihood ratio Lp: 
cos 9 m iss , the invariant mass of the ir° candidate, the 
p + candidate momentum, and the invariant mass of 
the 7T + 7r° pair used to make the p + candidate. The 
PDFs used in the likelihood ratio for the signal and 
background are determined from signal and B + B~ MC 
samples, respectively. 

The most powerful discriminating variable is B extra , 
defined as the sum of the energies of the neutral clus- 
ters not associated with the tag-B or with the signal tt° 
from the r + — > p + v mode, and passing a minimum en- 
ergy requirement (60 MeV). Signal events tend to peak 
at low -E ox tra- Background events, which contain addi- 



tional sources of neutral clusters, tend to be distributed 
at higher values. The signal region in data is kept blind 
until the end of the analysis chain when we extract the 
signal yield, meaning that we do not use events in data 
with Boxtra < 400 MeV during the selection optimization 
procedure and for the evaluation of background shapes. 

We optimize the selection requirements, including 
those on the purity V of the tag-B and the minimum 
energy of the neutral clusters, minimizing the expected 
uncertainty in the branching fraction fit. In order to es- 
timate the uncertainty, which includes the statistical and 
the dominant systematic sources, we run 1000 MC simu- 
lated pseudo experiments extracted from the background 
and signal expected B extra distributions for a set of pos- 
sible selection requirements, assuming a signal branching 
fraction of 1.8 x 10" 4 @|. 

Table [Til summarizes the signal selection requirements 
and Fig.[5]shows the B extra distribution with all the selec- 
tion requirements applied. The background events popu- 
lating the low B ext ra region are mostly semileptonic B de- 
cays for the leptonic modes. For the r + — s- it + v mode the 
background is composed mostly of charmless hadronic B 
decays and semileptonic B decays with a muon in the 
final state. For the r + —> p + v mode the backgrounds are 
charmed hadronic B decays, semileptonic B decays with 
a muon in the final state and a small fraction with a r. 

TABLE II: Optimized signal selection criteria for each r 
mode. 



Variable 


e+ 




7T + 




V 




> 10% 




Cluster energy ( MeV) 




> 


60 




TZ2 


< 0.57 


< 0.56 


< 0.56 


< 0.51 


COS 0TB 


< 0.95 


< 0.90 


< 0.65 


< 0.8 


Lp 






> 0.30 


> 0.45 



We use an extended unbinned maximum likelihood fit 
to the measured B cx tra distribution to extract the B + — > 
t + v branching fraction. The likelihood function for the 
N k candidates reconstructed in one of the four r decay 
modes k is 

e -(n s ,k+n b , k ) ( l 
Ck = II { n ^n(Ei, k ) + n b , k Vl(Ei,k)\, 

. (3) 

where n s k is the signal yield, rif, k is the background yield, 
E^k is the Bextra value of the i th event, V k is the PDF of 
signal events, and V k is the PDF of background events. 
The background yields in each decay mode are permitted 
to float independently of each other in the fit, while the 
signal yields are constrained to a single branching ratio 
via the relation: 

n s ,k = N B g x e k x B (4) 

where e k is the reconstruction efficiency of a particular 
r decay mode, and B is the B + — > t + v branching frac- 
tion. The parameters N g -s and are fixed in the fit 
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FIG. 2: -Ecxtra distribution in data (points with error bars) 
with all selection requirements applied and fit results overlaid. 
The hatched histogram is the background and the dashed 
component is the best-fit signal excess distribution. Plot (a) 
shows all r decay modes fitted simultaneously. Lower plots 
show the projection of the simultaneous fit result on the four 
analyzed r decay modes: (b) r + -» e + vv , (c) t + — > [i^vv , 
(d) r+ -> -k+u , (e) r+ -> . 



while S is allowed to vary. The reconstruction efficiencies 
efe, which include the r branching fractions, are obtained 
from MC-simulated signal events (see Table IIIip . Since 
the tag-B reconstruction efheiency is included in efc and 
is estimated from the signal MC, we apply a correction 
factor of i?data/MC = 0.926 ± 0.010 to take into account 



data/MC differences. This is derived from the ratio of 
the peaking component of the toes distribution for the 
hadronic tag-i? in data and in MC simulated events. 

The signal PDF is obtained from a high statistics sig- 
nal sample of MC simulated data. We use a sample of 
fully reconstructed events to correct the signal PDF for 
data/MC disagreement In addition to the reconstructed 
tag-S, a second B is reconstructed in the hadronic or the 
semileptonic decay mode using tracks and neutral clus- 
ters not assigned to the tag-B. In order to estimate the 
correction to the signal PDF, we compare the distribu- 
tion of .Ecxtra in this double tagged event sample from 
experimental data and MC simulations. The MC distri- 
butions are normalized to the experimental data and the 
comparison is shown in Fig. [31 We extract the correction 
function by taking the ratio of the two distributions and 
fitting it with a second order polynomial. 




Hadronic DT, E „ [GeV] 

' extra L J 



Hybrid DT, E [GeV] 



FIG. 3: -Ecxtra distribution for double tagged events. The 
"signal" B is reconstructed in hadronic decays (left plot) or 
semileptonic decays (right plot). Points are data and boxes 
are MC simulation. 



We determine the PDF of the combinatorial back- 
ground from the ibes sideband. The normalization of 
this component in the signal region is obtained by fit- 
ting the toes distribution after the selection has been 
applied. The shape of the peaking background is taken 
from B + B~ MC. The two background components are 
added together into a single background PDF. We esti- 
mate the branching fraction by minimizing — ln£, where 
C = Ii\ =1 Ck, and Ck is given in Eq. [3] The projections 
of the fit results are shown in Fig. [2j 

We observe an excess of events with respect to the ex- 
pected background level and measure a branching frac- 
tion of B(B+ -> t+v) = (1.83±g:||) x 1(T 4 , where the 
uncertainty is statistical. Table Mil summarizes the re- 
sults from the fit. We evaluate the significance of the 
observed signal, including only statistical uncertainty, as 
S = \j2\\i(C s j r b/ Cb), where £ s +b and Cb denote the 
obtained maximum likelihood values in the signal and 
background, and the background only hypotheses, re- 
spectively. We find S = 4.2a. 

Additive systematic uncertainties are due to the un- 
certainties in the signal and background -Ecxtra PDF 
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TABLE III: Reconstruction efficiency e, measured branching 
fractions, and statistical uncertainty obtained from the fit 
with all the modes separately and constrained to the same 
branching fraction. The r decay mode branching fractions 
are included in the efficiencies. 

Decay Mode e fc (xl(T 4 ) Signal yield ff(xl(T 4 ) 



\j + vv 



p+v 



4.1 ±9.1 0.35±o73 
1.12+ ' 90 



2.47 ± 0.14 
2.45 ±0.14 

0.98 ±0.14 17.1 ±6.2 3.691^ 

j+1.65 



12.9 ±9.7 
17.1 ± 6.2 
1.35 ±0.11 24.0 ±10.0 3.78 



-0.78 



■1.45 



combined 



62.1 ± 17.3 



1 oq+0.53 
-L.OO_q 49 



In summary, we have measured the branching frac- 
tion of the decay B+ — > t+v using a tagging algorithm 
based on the reconstruction of hadronic B decays us- 
ing a data sample of 467.8 x 10 6 BB pairs collected 
with the BA&4R detector at the PEP-II B-Factory. We 
measure the branching fraction to be B(B + — > t+v) = 
(1.83+^(stat.) ±0.24(syst.)) x 10~ 4 , excluding the null 
hypothesis by 3.8<r. (including systematic uncertainty). 
This result supersedes our previous result using the same 
technique 0. Combining this result with the other 
BABAR measurement of B(B + — > t+v) derived from a 
statistically independent sample [9(, we obtain B(B + — » 
t+v) = (1.79 ± 0.48) x 10~ 4 , where both statistical and 
systematic uncertainties are combined in quadrature. 



TABLE IV: Contributions to systematic uncertainty on the 
branching fraction. 



Source of systematics B uncertainty (%) 


Additive 


Background PDF 


10 


Signal PDF 


2.6 


Multiplicative 


Tag-B efficiency 


5.0 


B counting 


1.1 


Electron identification 


2.6 


Muon identification 


4.7 


Kaon identification 


0.4 


Tracking 


0.5 


MC statistics 


0.6 


Total 


13 



shapes used in the fit. To estimate the systematic uncer- 
tainty in the background PDF shape we repeat the fit of 
the branching fraction with 1000 variations of the back- 
ground PDFs, varying each bin content within its statis- 
tical uncertainty. We use the range of fitted branching 
fractions covering 68% of the distribution as systematic 
uncertainty yielding an overall contribution of 10%. We 
correct the systematic effects of disagreements between 
data and MC .©extra distributions for signal events us- 
ing a sample of completely reconstructed events in data 
and MC, as already described. To estimate the related 
systematic uncertainties, we vary the parameters of the 
second-order polynomial defining the correction within 
their uncertainty and repeat the fit to the B + — > t+v 
branching fraction. We observe a 2.6% variation that we 
take as the systematic uncertainty on the signal shape. 
Including the effects of additive systematic uncertainties, 
the significance of the result is evaluated as 3.8<r. 

Multiplicative systematic uncertainties on the effi- 
ciency stem from the uncertainty in the tag-i3 efficiency 
correction (5.0%), electron identification (2.6%), muon 
identification (4.7%), charged kaon veto (0.4%), and the 
finite signal MC statistics (0.8%). Table HVl summarizes 
the systematic uncertainties. The total systematic uncer- 
tainty is obtained by combining all sources in quadrature. 




0.2 0.3 
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500 1000 
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500 1000 
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FIG. 4: Top plot: Comparison between the measured 
B(B+ — > t+v) branching fraction (horizontal band) with the 
prediction of the 2HDM as a function of tan/3/m ff +, using 
exclusive (red/light gray) or inclusive (blue/dark gray) | V^t. | 
measurement. Bottom plots: 90% and 99% C.L. exclusion re- 
gions in the (m H + , tan ft) plane using the exclusive (left) and 
inclusive (right) measurements of \V u b\- 



Our measurement of the branching fraction B(B+ — > 
t+v) exceeds the prediction of the SM determined using 
the values of \V u b\ extracted from exclusive semileptonic 
events and from inclusive semileptonic events by 2.4er 
and 1.6cr, respectively. We also determine, separately for 
the exclusive and inclusive \V u b\ BABAR measurements, 
90% C.L. exclusion regions in the parameter space of 
the 2HDM - type II (to#+ , tan/3), where m H + is the 
charged Higgs mass and tan /3 is the ratio of the vacuum 
expectation values of the two Higgs doublets. We find 
that, taking \V u b\ from the exclusive measurement, most 
of the parameters space is excluded at 90% C.L. Using 
the higher value of \V u b\ from the inclusive measurement, 
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the constraints are less stringent but already set a lower 
limits at the TeV scale for high tan/3. The same implica- 
tions on 2HDM are supported by a recent BABAR study 
of the B(B -> D^tv) decays QJ]. Fig. Ushows a com- 
parison between the measured B(B + — > t + v) branching 
fraction with the prediction of the 2HDM as a function of 
tan/3/m#+ and the exclusion plots in the (mn+ , tan 0) 
plane for the exclusive and inclusive measurements of 
\V ub \ . 
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